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Lewis Research Center 

SUMMARY 

A divergent-flow contact-ionization electrostatic thrustor is described herein. This 
thrustor performed well at thrust levels between 0. 89 and 1. 56 milline wtons, a range 
that is applicable to satellite attitude -control and station-keeping missions. Power effi- 
ciencies up to 45 percent (excluding vaporizer and neutralizer powers) were obtained for 
operation at specific impulses from 5000 to 8000 seconds. No damage due to primary 
ion impingement was detected in about 200 hours of running. Accelerator drain currents 
were typically less than 2 percent of the beam current. Computer studies indicated an 
accelerator electrode lifetime in excess of 20 000 hours. 

Experimental data verified the computer predictions for electrostatic beam deflec- 
tion up to about 15°. This deflection produced a component of thrust normal to the 
center plane of up to 30 percent of the total thrust of an undeflected beam. 


INTRODUCTION 

Contact-ionization thrustors using cylindrical electrode geometries are investigated 
analytically in reference 1. Most experimental effort to date has been directed toward 
converging flow and plane parallel flow. This report presents theoretical and experi- 
mental evaluation of a thrustor based on diverging flow between coaxial cylinders 


Part of the material, " Experimental Performance of a Low-Thrust, Divergent- 
Flow, Contact-Ionization Electrostatic Thrustor" by John F. Staggs and Walter C. 
Lathem, was presented at the AIAA Second Propulsion Joint Specialist Conference, 
Colorado Springs, Colo. , June 13-17, 1966, and was published in the Journal of Space- 
craft and Rockets, vol. 4, no. 5, May 1967, pp. 610-615. 
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(refs. 1 and 2), as shown in figure 1. The thrustor develops a thrust within a range suit- 
able for attitude control and station keeping of satellites in the 200- to 700-kilogram mass 
range (ref. 3). With electrostatic beam deflection, a thrustor can perform both 
attitude -control and station-keeping functions without thrustor gimbaling. 

Two important thrustor performance requirements are long operational life and high 
efficiency. For satellite missions of several years duration, long life is more impor- 
tant, since these thrustors will require only a small portion of the satellite power supply. 

To obtain long operational life, several components must be considered. The life- 
times of certain components, such as the porous tungsten ionizer and the ionizer heater 
element, are dependent to a large degree on the state of the art. Accelerator grid life, 
however, is strongly dependent on the optical design of the thrustor; that is, the optics 
will determine the rate of accelerator erosion due to charge -exchange ion impingement. 
The primary advantage of a diverging ion beam is that both the ion density and the neutral 
atom density are reduced in the accelerator aperture area, and the charge exchange ion 
formation rate is thus reduced. (The charge -exchange interaction rate is directly pro- 
portional to each of these quantities. ) 

In this report, the design, fabrication, and testing of a divergent-flow contact- 
ionization thrustor is described in detail. 

THEORETICAL CONSIDERATIONS 
Electrode Geometry 

The electrode system used in the experimental thrustor was designed by using an 
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Figure 2. - Cross section of ionizer, electrodes, and optical characteristics from digital computer 
program. 
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Figure 3. - Ideal potential distribution. 

electrolytic tank analog (ref. 4) and a digital computer program (ref. 5). The goal was 
to provide a reasonably high current density with no primary ion impingement on the ac- 
celerator. A secondary goal was simplicity of electrode fabrication and alinement. A 
cross-sectional sketch of the electrode shapes arrived at is shown in figure 2. Also 
shown are typical ion trajectories and equipotential lines for the Poisson solution ob- 
tained with the digital program. The ideal potential distribution is shown in figure 3. 

The computed variation of current density with total voltage can be compared with 
basic theory by considering the theoretical equation for space-charge limited flow be- 
tween coaxial cylinders given in reference 2 as 
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(All symbols are defined in appendix A. ) When cesium is used, this equation becomes 
(in SI units) 


h 


4.74X10" 9 



A/sq m 


In the experimental configuration r o = 3. 0 millimeters and r = 6. 0 millimeters 
(/3 = 0. 279). The solution of the preceding equation using these values yields the upper 
curve of figure 4. Digital solutions for the final electrode configuration are also shown 
for two values of the focuser-accelerator spacing b. The ionizer-focuser spacing a 
was held constant at 0. 127 millimeter, and the accelerator aperture was maintained at 
3. 0 millimeters. The computed current density for the thrustor configuration is lower 
than the ideal theoretical values for two reasons. First, the overall electric field is re- 
duced by the aperture effect. Second, the presence of a focusing electrode acts to reduce 
the electric field at the ionizer. This latter effect is apparent from inspection of 
figure 5 which shows the results of varying the spacing a, from 0. 0 to 0. 191 millimeter 



Figure 4. - Comparison of theoretical 
and digital computer program current 
densities as functions of total accel- 
erating voltage for two values of accel- 
erator focuser spacing. 



Figure 5. - Effect of ionizer-focuser spacing 
on emitter current density from computer 
program. Focuser-accelerator spacing, 

2 . 14 millimeters. 
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Figure 6. - Fraction of accelerator aperture occupied 
by ion beam as function of ionizer-focuser spacing. 

while holding b constant. 

The fraction of accelerator aperture occupied by the beam as a function of the 
ionizer-focuser spacing a is shown in figure 6. Figure 6 shows that beam divergence 
is a sensitive function of this spacing. 


Accelerator Electrode Life 

The useful life of an ion thrustor accelerator electrode is limited by the sputtering 
rate of the electrode material due to impinging ions. The impinging ions can be of two 
types, primary beam ions not properly focused or ions formed in the electrode aperture 
area by charge exchange of primary ions and neutrals from the emitter. 

Figure 2 (p. 3) shows that, theoretically, the divergent flow thrustor beam occupies 
only about 60 percent of the accelerator aperture. This precludes the possibility of ero- 
sion due to primary ions. There will, however, be erosion from charge -exchange ions. 
As stated earlier, one of the theoretical advantages of the divergent-flow concept is the 
reduced ion and neutral density in the aperture area, which should result in a reduction 
in the rate of charge -exchange erosion, thereby extending thruster lifetime. In order to 
predict the accelerator charge-exchange erosion and thrustor lifetime, an analysis was 
made utilizing both the electrolytic tank analog and the computer programs. This anal- 
ysis yielded the calculated erosion patterns shown in figure 7. The parameters used for 
the study are given in figure 7, which is drawn for 20 000 hours of operation at these 
conditions. For details of the calculation and further discussion of the techniques used 
in this study, see appendix B and also references 4 and 6. 
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Figure 7. - Theoretical accelerator erosion pattern due to charge-exchange ion impingement 
after 20 000 hours at 163 amperes per square meter. Total beam current, 16. 3 milliamperes; 
neutral atom loss, 1.0 percent of total beam current. 


EXPERIMENTAL APPARATUS 
Thrustor Design 

A cross-sectional sketch of the experimental thrustor configuration is shown in 
figure 8 and photographs are shown in figure 9. The overall size of the thrustor (in- 
cluding a grounded shielding screen not shown in the photographs) is 2. 5 by 7. 5 
by 14. 0 centimeters. The total weight of the thrustor, excluding the vaporizer, is ap- 
proximately 0. 15 kilogram. 

The accelerator electrode was made of copper, so that back -sputtered material 


— i— r — i 


4. 0 cm 



^-Accelerator electrode (copper) 
✓-Focus electrode (molybdenum) 


,-Boron 

nitride 


I fjfll -nca ici •• y\ \ 

; r • Stri P 


^^JfiiFibrous thermal ^,7 
fepSS! insulation 


Stainless-steel mounj plate 
1 75 cm 


6 


Figure 8. - Cross section of thrustor showing overall dimensions of support- 
ing structure and materials used. 



Focus er electrode — ^ 



065-273 


(a) Partly disassembled. 



067-1796 

(b) Assembled. 

Figure 9. - Thrustor. 

would not clog the porous ionizer. It was also made in two halves to allow electrostatic 
beam deflection by varying the potential on each half. The ionizer, electrodes, and 
heater lead-in rods were all mounted from the stainless-steel mount plate (fig. 9(a)), 
which was maintained at ionizer and focuser electrode potential. The focuser was made of 
molybdenum. Boron nitride sideplates confined the fibrous thermal insulation material. 
In operation, a stainless -steel grounded screen surrounded the thrustor to prevent stray 
electrons from reaching the thrustor. 

The ionizer assembly is shown in figure 10. The ionizer manifold consists of a U- 
shaped channel of 0. 127 -millimeter -thick tungsten sheet electron-beam-welded to the po- 
rous tungsten ionizer. End plates, feed tube, and support tabs are of tantalum. Another 
U-shaped channel of tungsten is welded to the back of the manifold to form an enclosure 
for the heater strip (fig. 8). Heating the ionizer is accomplished by radiation from the 
strip. The ionizer assembly, which is approximately square in cross section (with each 
side 2. 2 mm long), has an ion emitting surface area of 1. 0 square centimeter. This 
convex emitting surface has a radius of 3. 0 millimeters, a chord length of 1. 52 milli- 
meters, and an overall length of 6. 35 centimeters. 
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(b) Exploded. 

Figure 10. - Ionizer assembly. 

Test Facility 


A schematic drawing of the test setup is shown in figure 11. The main chamber of 
the vacuum facility is 4. 5 meters long by 1. 5 meters in diameter and contains a liquid- 
nitrogen- cooled condenser. The facility was typically pumped to a pressure of from 
10 -D to 10 torr (approx. 10 to 10"° N/sq m) during tests by using three 0. 9-meter 
oil diffusion pumps. Actual pressure depended on operating conditions. A 0. 9-meter- 
diameter gate valve separated the main chamber from the thrustor chamber, which was 
approximately 1 meter long by 1 meter in diameter. A more detailed description of the 
vacuum facility is given in reference 7. All thrustor electrical connections were made 
through ceramic insulators at the rear of the thrustor chamber. 

An isolation screen in the exit plane of the thrustor was maintained at electrical 
ground to prevent electrons from reaching the electrical leads in the thrustor chamber. 


Instrumentation 

The meters used in measuring beam currents and accelerator drain currents had an 
accuracy of ±1. 0 percent of full scale. They were provided with multiple shunt circuits 
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to allow full-scale current readings from 0. 1 to 30 milliarnperes. 

A molybdenum -tipped probe assembly was used to investigate beam dispersion, to 
obtain a check of the total ion beam current obtained from the meter readings, and to 
check the amount of beam vectoring when a differential voltage was applied to the two ac- 
celerator halves. The assembly of 15 buttons, each 12. 7 millimeters in diameter, was 
located 35. 6 centimeters downstream of the thrustor. As the assembly traversed the ion 
beam, the currents produced by ions striking the probe buttons created a voltage drop 
across 100 -ohm terminating resistors. A system capable of handling 20 data points per 
second was used to sample the 15 buttons plus a position indicator signal in rotation. In 
this way, each button was sampled approximately 75 times per minute. The probe 
traveled about 40 centimeters per minute to give approximately two data points per centi- 
meter of travel for each button. A more detailed discussion of molybdenum-tipped probe 
operation is given in reference 8. 
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RESULTS AND DISCUSSION 


Ionizer Current Density and Critical Temperature 

During testing of the thrustor, the ionizer-to-focuser spacing was held at 0. 127 mil- 
limeter, and two values of focuser-to-accelerator spacing were used. A comparison of 
the experimental current densities with computer predicted values at various cesium 
flow rates is given in figure 12 for b = 1. 52 and 2. 14 millimeters. The maximum cur- 
rent density obtained experimentally was 160 amperes per square meter at a field 
strength of 3. 0x10 volts per meter. As predicted by the computer program, varying 
from 1. 5 to 3. 0 with constant <& did not significantly change the experimental 
current density. 

A high-quality porous tungsten designated M E-4 M in reference 9 was used throughout 
the program. The samples were fabricated from 1- to 5 -micron-diameter spherical 
powder and were about 79 percent of theoretical density. Critical temperatures for vari- 
ous ionizer current densities were obtained from pellet samples and are reported in ref- 
erence 10. These values are compared with results of this program in figure 13. Also 
included in the figure is the zero-field Taylor- Langmuir curve for solid tungsten taken 
from reference 11. In the pellet evaluations of reference 10, the critical temperature 



(a) Focuser-accelerator spacing, (b) Focuser-acceleratnr spacing, 
1.52 millimeters. 2. 14 millimeters. 

Figure 12. - Comparison of experimental current densities with computer 
predicted values, lonizer-focuser spacing, 0. 127 millimeter. 



Figure 13. - Comparison of critical 
temperature measurements with 
those of reference 10 and with the 
zero-field Taylor-Langmuir curve 
for solid tungsten (ref. 11). 
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points were taken as points of minimum neutral fraction, whereas the data reported 
herein reflect temperatures at which the current density was 95 percent of its maximum 
value. At these temperatures, the accelerator drain currents were near a minimum. 
Also, since the ionizer temperature was somewhat nonuniform along its length, the crit- 
ical temperature was taken to be the coolest point on the ionizer surface. 


Accelerator Drain Currents 


Typical data of the ratio of total accelerator drain current to beam current Jg 
as a function of total accelerating voltage $ are presented in figures 14 and 15. In fig- 
ure 14, data are compared for three different flow rates. A possible explanation of the 
reversal of the slope of the curves as the voltage is increased might be field enhance- 
ment of thermionic emission from the cesium-coated accelerator (ref. 12). In figure 15, 
results obtained with a background pressure of 7. OxlO - ^ torr (approx. 9x10"^ N/sq m) 

are compared with those obtained at the same mass -flow rate but with a background 
-5 -3 

pressure of 6.0x10 torr (approx. 8x10 N/sq m) of oxygen. When oxygen is admitted 
to the system, it may change the thermionic emission characteristics of the cesiated 
copper accelerator surface as well as have a decarburizing or oxidizing effect on the 
ionizer (ref. 13). In every instance the addition of oxygen reduced the accelerator drain 
currents. The exact mechanism of the change is not fully understood, but it is believed 



Figure 14. - Ratio of accelerator drain current to 
beam current as function of total accelerating 
voltage for several mass-flow rates. Background 
pressure, 8xl0~7 torr (approx. HT^N/sq m). 
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Figure 15. - Effect of oxygen 
atmosphere on accelerator 
drain currents. Calcu- 
lated mass-flow rate, 

8. 9x10 kilogram per 
second. 
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that the reduction in drain currents occurs too rapidly after the addition of oxygen (on 
the order of 1 min) to be caused by decarburization of the ionizer and is probably due to 
a combination of the other effects. 

Possible sources for accelerator drain currents include direct impingement, 
charge-exchange ion impingement, leakage currents across cesium-coated insulators, 
thermionic electron emission from a hot cesium-coated accelerator, and secondary 
emission resulting from ion impingement. 

The digital computer program predicted no direct ion impingement. After a total of 
200 hours of thrustor operation at an average beam current of 6. 2 milliamperes, the 
accelerator grid showed no observable erosion in the areas which might be struck by 
ions from the ionizer. 

Some charge exchange, however, is inevitable. The magnitude depends not only on 

the neutral fraction but also to a large extent on the optics of the system. A diverging 

system reduces the ion density and the neutral atom density in the accelerator aperture 

region and thus minimizes the rate of charge -exchange interactions. From the charge- 

exchange study mentioned previously, the total charge -exchange impingement rate was 

12 

calculated to be about 10 ions per second, or approximately 0. 001 percent of a 16. 0- 
milliampere ionizer current. Charge -exchange, then, is probably a negligible contri- 
butor to the drain currents observed. 

Since ion impingement (both direct and from charge exchange) appears to be negli- 
gible, secondary electron emission would also be very small. All insulators were to 
the rear of the thrustor and were therefore not subject to cesium coating. Leakage cur- 
rents across these insulators were negligible. Based on these arguments, it is likely 
that nearly all the experimental drain currents can be attributed to thermionic electron 
emission from the upstream face of the accelerator, which probably became coated with 
cesium atoms from the ionizer. 

Two important factors affecting the neutral fraction besides the characteristics of 
the tungsten are (1) temperature uniformity over the* ionizer and (2) matching the arrival 
rate profile of cesium of the emitting surface with the ion extraction rate profile from 

< 1.52 mm ► 

r-Downstream surface radius 
\ of_curvature, 3. 0 mm 


r -Pores on these surfaces r-Upstream surface radius 

'^closed by electron-beam welder^ __ ^,of curvature, 1. 27 mm 
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Figure 16. - Cross section of porous tungsten ionizer showing surfaces that have been 
closed with electron -beam welder. 


that surface. 

The first factor can be handled by designing the ionizer and heater to assure tem- 
perature uniformity and/or operation with the coolest portion of the ionizer above the 
critical temperature. 

The second factor is taken care of by contouring the back (upstream face) of the po- 
rous tungsten. The final shape (fig. 16) was obtained by setting up an analog of the flow 
through the porous tungsten by using the electrolytic tank. In this figure, the dashed 
lines show the contour of the porous surfaces and the solid lines represent surfaces 
where the pores have been ” washed" closed with an electron-beam welder (short dura- 
tion heating, so that a thin surface layer is melted) . 


Optical Characteristics and Thrust Vectoring 


Optical characteristics of the thrustor were investigated by beam probe measure- 
ments. Figure 17 shows contour plots of the current -density variation through a cross 
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Molybdenum probe assembly 
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(b) Deflection voltage of 1. 5 kilovolts applied to split accelerator. Beam deflected about 15°. 

Figure 17. - Current contour map from molybdenum button probe measurements showing lines of equal current density. Total current calcu- 
lated from this figure, 14. 7 milliamperes. 


13 



section of the beam 35.6 centimeters downstream of the ionizer. Figure 17(a) shows an 
undeflected beam, while figure 17(b) shows the effect of a 1. 5-kilovolt deflection voltage. 
This deflection voltage was obtained by increasing the voltage on one accelerator half by 
750 volts and decreasing the voltage by the same amount on the other half. 

While these data were being taken, the metered beam current was 9. 4 milliamperes. 
A rough integration of this contour map yielded a current of about 14. 7 milliamperes. 

The difference is due both to secondary emission from the probes and to inaccuracies in 
drawing the map. Although these current values do not show good agreement, the quali- 
tative profile and magnitude of beam deflection obtained agreed with theoretical predic- 
tions. For example, it was determined from the contour map that 90 percent of the 
beam lies within a dispersion angle 9 of about 15°. This value agrees well with the 
18° angle predicted by the digital computer solution (fig. 2, p. 3). The double peaked 
profile of the beam had also been predicted from the computer results. 




Figure 18. - Effects of applying differential voltage A (p to two halves of split 
accelerator. Solutions obtained from digital computer; A <p = <pi~ (fty 
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The beam deflection angle a was analyzed theoretically, and typical results are 
shown in figure 18 for a differential voltage A <p of 1. 5 and 0 kilovolt. The deflection 
angle or (fig. 19) is approximately 1. 0° for every 100 volts of deflection voltage. (This 
angle is not to be confused with the dispersion angle 6 , which is a measure of the nor- 
mal spread of the beam.) Experimental results for beam deflection are compared with 
theory in figure 20. The deflection angles determined from the probe contour maps are 
shown as data points and are in good agreement with theory. 

The calculated thrust component normal to the center plane T as a fraction of T 

y ^ 

for a = 0 is plotted against the calculated deflection angle in figure 21 (as in fig. 19, 

T x = T reg for a = 0) . With a deflection of 15°, Ty is about 30 percent of (T x ) a _Q 



Center 

plane 


Figure 19. - Cross section of thrustor showing deflection angle a and 
thrust vectors. 



Figure 20. - Beam deflection angle as function of 
deflection voltage. 



thrust vector, a 

Figure 21. - Thrust component normal to center 
plane (as fraction of undeflected total thrust) 
as function of deflection angle. 
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Overall Performance 


The thrustor yielded overall performance which approximates the requirements for 
certain satellite attitude -control and station-keeping missions in which a thrust range of 
0. 445 to 1. 78 millinewtons (0. 1 to 0. 4 mlb) is needed (ref. 3). Stationary-satellite 
missions require control systems that are lightweight, efficient, and long lived. In the 
continuous corrections mode of operation, approximately 0. 445 millinewton (0. 1 mlb) of 
thrust is needed for every 160 kilograms of satellite weight (ref. 3). Figure 22 shows the 
calculated thrust as a function of beam current for three different specific impulses 
(thrust calculated from measured beam currents and net accelerating voltages) . This 
variation in specific impulse was achieved by changing $/<P-g while maintaining con- 
stant 4>. 

Figure 23 shows power efficiencies as a function of specific impulse, calculated 
from 




P B + P h + P A 


These efficiency values do not include vaporizer or neutralizer powers because these 
components were not optimized for the tests. Accelerator drain power losses are in- 
cluded. Data are presented for accelerator-to-focuser spacings of 1. 52 and 2. 14 milli- 
meters. 



Figure 22. - Calculated thrust as 
function of experimental beam 
current for three values of 
specific impulse. 
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Figure 23. - Thrustor power efficiency as function 
of specific impulse. 
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Figure 24. - Power-thrust ratio as function of spe- 
cific impulse. 


Figure 24 shows the variation of power-thrust ratios with specific impulse at dif- 
ferent beam current levels. The dashed curve represents a 100-percent-efficient thrus- 
tor. 

Figures 23 and 24 show that efficiency could be improved considerably. However, 
as previously noted, the more important factor for satellite applications is lifetime. 

Thus far, over 200 hours of testing have indicated no problems with accelerator erosion. 
More testing would be required to verify the accelerator erosion predictions fully. Op- 
timization of components such as the ionizer heater assembly, electrode support struc- 
tures, and thermal insulation package should improve overall performance. 


CONCLUDING REMARKS 

A contact-ionization electrostatic thrustor was designed, fabricated, and tested. 
Thrust levels between 0. 89 and 1. 56 millinewtons (0. 2 and 0. 35 mlb) were produced at 
power efficiencies up to 45 percent in the specific-impulse range from 5000 to 8000 sec- 
onds. The optical characteristics were such that primary ion impingement on the accel- 
erator electrode was negligable or nonexistent, and drain currents were about 2 percent 
or less of the beam current. Computer studies and numerical analysis predicted elec- 


17 


trode lifetimes in excess of 20 000 hours. 

Electrostatic ion beam deflection worked well for producing a thrust component nor- 
mal to the center plane. By deflecting the ion beam 15°, a thrust is produced in a direc- 
tion perpendicular to the thrustor center plane which is about 30 percent of the original 
total thrust. With this technique, both attitude -control and station-keeping functions 
could be performed with the same thrustor, and the need for mechanical gimbaling would 
thus be eliminated. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 16, 1967, 

120-26-02-06-22. 
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APPENDIX A 


SYMBOLS 


a 

ionizer-focuser spacing 

w 

width of emitter chord 

b 

focuser -accelerator spacing 

X 

normal distance from emitter 

f 

ratio of neutral atom emission 
rate to ion emission rate 


surface to center of charge - 
exchange volume 


from emitter 

Y 

normal distance from plane of 

I 

sp 

specific impulse 


symmetry of thrustor to 
center of charge -exchange 

J 

current 


volume 

j 

current density 

a 

beam deflection angle 

k 

Boltzmann constant 

(3 

function of r/r Q (ref. 2) 

m 

mass of atom 

1 ?P 

power efficiency 

m 

mass -flow rate 

6 

beam dispersion angle 

N 

charge -exchange ion formation 

M 

arrival rate at charge- 


rate 


exchange volume 

P 

power 

V 

emission rate from ionizer 

Q 

cross section for charge - 

P 

number density 


exchange interaction 

$ 

total accelerating voltage 

q 

charge on electron 

<P 

voltage 

r 

radius of outer cylinder of 
coaxial cylinder electrode 
system 

Subscripts: 

A 

accelerator 

r o 

radius of inner cylinder (emit- 

B 

beam 


ter) of coaxial cylinder 

E 

emitter 


electrode system 

h 

ionizer heater 

T 

temperature 

j 

j plane parallel to yz -plane 

^res’ 1 
T x’ T yJ 

thrust vectors 

1U 

k k in plane parallel to xz -plane 

Superscripts: 

V 

charge -exchange volume 

0 

neutral atom 

V 

mean velocity 

+ 

positive ion 
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APPENDIX B 


CHARGE -EXCHANGE EROSION CALCULATIONS 


Calculating the effect of charge -exchange ion erosion on an ion thrustor accelerator 
requires knowledge of 

(1) The rate charge -exchange ions are formed in the primary ion beam region 

(2) The arrival rate of charge-exchange ions as a function of energy and angle of in- 

cidence for each element of accelerator electrode area 

(3) The sputtering rate of the accelerator material 

The configuration used for this figure is shown in figure 25. It will be helpful to refer to 
this figure in following the analysis to be presented. 



Figure 25. - Configuration used for charge-exchange erosion analysis. Charge-exchange ion trajectories started at centers of mesh squares. 


The interaction cross-section equation is 





(Bl) 


The numbers j and k identify a volume element in the beam region. 

The interaction cross section Q for cesium ions and cesium atoms is given by 
Marino (ref. 14) as 


Q, , = (26. 8 - 1. 46 In <p. J 2 xi0' 20 

J? & J? ^ 


(B2) 
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The voltage <p. , is the difference between emitter potential and the potential at the 

center of the j, k charge-exchange volume and is determined from the digital computer 

solutions. The ion arrival rate pt , is also obtained from the digital computer solu- 

1 > K 

tion. 

Obtaining the neutral atom density p? . requires a neutral atom emission rate 
0 1> K 

v from the emitter, which can be expressed as a fraction f of the ion emission rate 
v + ~, that is, 


-2 +0 

A value for f of 10 is assumed. The quantities v and v are assumed constant 
across the emitter. The atom arrival rate is then calculated from the equation (ref. 15) 



where 


v° 



The arrival rate of impinging charge -exchange ions on the accelerator is determined 
by tracing out trajectories on the electrolytic tank analog. The Poisson potential distri- 
bution from the computer solution is established in the electrolytic tank. The exhaust 
beam is divided into a number of mesh squares (fig. 25, p. 20). Charge -exchange ion 
trajectories are then started with zero initial velocity from the centers of these mesh 
squares and traced on an xy-plotter. Figure 25 shows that the accelerator electrode has 
been divided into elemental segments on the downstream face. (No ions were found to 
impinge on the upstream face.) All ions formed in a particular mesh square are as- 
sumed to begin at the center of that square and to have the same energy when striking 
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the grid. Furthermore, it is assumed that they are equally distributed over the ele- 
mental grid segment and strike at the same angle of incidence. These assumptions are 
reasonable since the mesh squares and accelerator electrode segments are compara- 
tively small. With these assumptions, the next step is to determine the sputtering dam- 
age for each elemental segment. Sputtering yields as a function of energy and angle of 
incidence for cesium ions (amu, 132. 9) incident on copper were not available. Instead, 
data for xenon (amu, 131. 3) on copper from reference 16 was used. 

In summary, the charge -exchange ion formation rate from equation (Bl) and the re- 
sults of the electrolytic tank trajectory plots (fig. 25) yielded a solution for the distri- 
bution of ions striking the downstream face of the accelerator electrode. Then with the 
use of the sputtering data of reference 16 an erosion pattern was obtained. An example 
of an erosion pattern after 20 000 hours is shown in figure 7 (p. 6). 
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